Thermodynamics of the anisotropic two-channel Kondo problem 



O 

o 

(N 
o 

o 

o 



I 

o 
o 



o 

O 



o 



o 
o 



X 



Gergely Zarand^'^, Theo Costi'^, Andres Jerez^, and Natan Andrei''^ 
^ Lyman Physics Laboratory, Harvard University, Cambridge, MA 
^Research Group of the Hungarian Academy of Sciences, Institute of Physics, TU Budapest, H-1521 
^ Institut Laue-Langevin, 6 rue Jules Horowitz, B.P. 156, 38042 Grenoble Cedex 9, France 
Center of Materials Theory, Rutgers University, Piscataway, NJ 08854-0849, USA 

(February 6, 2008) 

We construct and solve numerically the thermodynamic Bethe Ansatz equations for the spin- 
anisotropic two-channel Kondo model in arbitrary external field h. At high temperatures the specific 
heat and the susceptibility show power law dependence. For h ^ and at temperatures below the 
Kondo temperature Tk a two-channel Kondo effect develops characterized by a Wilson ratio 8/3, 
and a logarithmic divergence of the susceptibility and the linear specific heat coefficient. Finite 
magnetic field, h > drives the system to a Fermi liquid fixed point with an unusual Wilson ratio 
which depends sensitively on h. 

PACS numbers; 75.20.Hr, 71.10.Hf, 71.27.-|-a, 72.15.Qm 



I. INTRODUCTION 

The two-channel Kondo model (2CIKM) attracted a lot 
of interest during the past few years.El This intense inter- 
est is mostly triggered by the rather unusual properties 
of this model: The existence of finite residual entropy in 
zero external field, the vanishing of single particle scat- 
tering amplitude at T = temperature, and the loga- 
rithmic singularity of various thermodynamic quantities. 
All these unusual properties appear due to the presence 
of a hidden and conserved flavor quantum number of the 
electrons. 

Several physical systems have been proposed to realize 
the 2CKM. An unambiguous realization of the 2CKM 
is provided by dilute. Uranium and Cerium-based heavy 
fermion compounds.^ In these alloys the combination of 
strong spin-orbit interaction and crystal field symmetry 
effects leads to an effective 2CKM. 

It has been also suggested that non-commutative tun- 
neling centers may form two-level systems (TLS) and re- 
alize the 2CKM.B In this case the localized spin is re- 
placed by the position of the tunneling center, while 
the angular momenta of the conduction electrons play 
the role of the electron spin in the 2CKM. The elec- 
tron spins in the TLS problem play the role of silent 
flavor indices. Though there have been several concerns 
raised concerning the microscopic structure of the TLS 
and the possibility of observing the two-channel Kondo 
behaviorJ3'El many experiments may be consistently in- 
terpreted in. terms of such dynamical two-channel Kondo 
impurities. ETlI 

Another realization of the 2CKM is provided W quan- 
tum dots near to their charge degeneracy point .u In this 
case the charging states of the dot replace the impurity 
spin states and they couple to the position variable of 
the conduction electrons. Again, the spin of the elec- 
trons acts as a flavor variable. Though it appears to 
be extremely difhcult to observe experimentally the non- 
Fermi liquid behawior associated with the two-channel 
Kondo fixed point some fingerprints of the two-channel 



Kondo effect have been observed iii-|the capacitance of 
such semiconducting quantum dots.Ej 

The above physical realizations of the 2CKM have the 
common feature that the 2CKM that describes them is in 
all cases generically strongly spin-anisotropic, and breaks 
the full SU(2) spin symmetry. While this SU(2) sym- 
metry breaking is known Xp be irrelevant at the two- 
channel Kondo fixed point jll it affects both qualitatively 
and quantitatively the properties of the model at ener- 
gies around and above the Kondo temperature, Tk, and 
has to be taken into account to make comparison with 
experiments. 

In the present paper we generalize the Bethe Ansatz 
results of Ref . [ to compute both analytically and nu- 
merically the thermodynamics of the anisotropic 2CKM 
at arbitrary temperatures and magnetic fields for vari- 
ous values of the anisotropy. Below the Kondo temper- 
ature, Tff, we find that the thermodynamic properties 
of the model are very similar to those of the isotropic 
2CKM. This observation is in full agreement with the ir- 
relevance of the anisotropy at the 2CKM. However, above 
Tk the behavior of the model depends crucially on spin 
anisotropy. We find that above Tk all thermodynamic 
quantities display a power law behavior with an exponent 
determined by the value of the anisotropy. Furthermore, 
we find, that at finite magnetic fields the system flows to 
an unusual Fermi liquid fixed point with an anisotropy- 
dependent Wilson ratio. 

The paper is organized as follows. In Section || we in- 
troduce the 2CKM. In Section III we construct the ther- 



modynamic Bethe Ansatz equations fo r th e anisotropic 
2CKM that we analyze in detail in Sec. IV, Some details 
of the computation are given in the Appendix. 



II. MODEL 

The anisotropic 2CKM consists of a spin 5o = 1/2 im- 
purity that couples dynamically to the spin of the conduc- 
tion electrons through a strongly anisotropic exchange 
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interaction, and is described by the foUpwing first quan- 
tized Hamiltonian {h = = fJ-B = ^)'-^ 

j — 1 a—x.y,z 

+ h[gS^,+g' Y^Sj) . (1) 
j 

In this equation Xj is the coordinate of the j'th conduc- 
tion electron, — N x f is the total number of electrons, 
and the Sj's (j = 0, .., A^e) denote the spin 1/2 operators 
of the electrons. Electrons also have a conserved flavor 
quantum number, m = {!,..,/}, implicit in Eq. (^. In 
the following we discuss the case / > 2 as well though 
for the physical realizations mentioned in the Introduc- 
tion/ ~ 2. (The f — I case has been analyzed in Refs. [ 
11,0.) 

The meaning of the various terms in this Hamiltonian 
depends on the particular physical realization. The ex- 
ternal field h couples to the local moment and the conduc- 
tion electron spins with different g- factors. In the TLS 
problem and the quantum dot case, the external field h 
represents the asymmetry of the TLS and the splitting 
between the two charging states of the quantum dot, re- 
spectively, and g' = 0. For heavy fermion systems, de- 
pending on the particular realization, h can describe an 
external magnetic- or strain field, and may also couple to 
J2j Sj- In the present paper we mostly study the local 
susceptibility, corresponding to g' = 0. 

In the original formulation of the Bethe Ansatz, sim- 
ilar to tlie, conformal field theory solution of the Kondo 
problcmjlj the impurity spin is 'fused' with the spin de- 
grees of freedom of the conduction electrons, and the ex- 
ternal field couples to the total spin, corresponding to 
g' = 1. Thus the Bethe Ansatz calculates the impu- 
rity contribution to the global susceptibility with g' = g. 
Lowenatein made a remarkable attempt to treat the g' 
1 case,Efl however, his results were not fully conclusive .cfl 
In the present paper we use a different strategy to cope 
with the g' — situation, and show that, similar to ttte 
case of the single channel anisotropic Kondo problem,L3 
the g' = local susceptibility (i.e. the susceptibility aris- 
ing when the field couples only to the impurity spin) is 
simply proportional to the global susceptibility. There- 
fore, after determining the appropriate normalization fac- 
tor, we are able to extract the exact local susceptibility 
from the Bethe Ansatz calculation with g — g' ■ 

In our analysis we shall assume that Jz Jx — Jy — 
J_L and thus the Hamiltonian has a U{1) spin symme- 
try. This is true for the quantum dot and heavy fermion 
realizatipns, however, for the TLS problem Jz ^ Jx > 
Jy — On Fortunately, it is not necessary to treat the gen- 
eral Bethe Ansatz equations in order to determine the 
universal features of the TLS Hamiltonian in the TLS 
case either: Under scaling the terms Jx and Jy, describ- 
ing electron-assisted tunneling processes, become rapidly 
equalB and therefore it suffices to consider an equivalent 



effective model with = = J_l = to deter- 

mine the thermodynamics. This effective model has U{1) 
orbital symmetry and can be much more easily analyzed 
by Bethe Ansatz techniques than the fully anisotropic 
one. 

In constructing the wave function we implicitly have 
to introduce a further electron-electron interaction 

^el-el^ E E S{x,-X,)J^SrSf . (2) 

l<i<j<Ns a=x.y,z 

Since electrons move with the same velocity and in the 
same direction, this interaction does not modify the ther- 
modynamics in a crucial way for Jz > J±. In particular, 
we have shown that the specific heat of the electrons is 
unaffected by this interaction. However, it does rescale 
the g-factor in Eq. (|l|): g ~* g°^ . Fortunately, we shall 
be able to compensate this effect by the renormalization 
of g-factors in the magnetic field term. 

For Jz < J± the artificial electron-electron interaction 
apparently generates a gap in the spin sector, and there- 
fore it is impossible to capture the Kondo effect there 
within the algebraic Bethe Ansatz approach. Therefore, 
unfortunately, we have to restrict our discussion to the 
case Jz > J±, and our calculations are not directly ap- 
plicable to the quantum dot case. 

In Eq. |l], only a term ^ hS^ has been included, which 
acts as a local field in the z direction. In reality, how- 
ever, the field can also point into the 'perpendicular' di- 
rection corresponding to a term ^ h^S^. In the TLS 
case this term describes spontaneous tunneling between 
the TLS positions. The general effect of the latter term 
is very similar to that of hS^, but quantitatively it be- 
haves somewhat differently. For a TLS, e.g., in a tKpical 
situation h is believed to be much larger than h^p and 
furthermore, its effective value, is orders of magnitude 
reduced by polaronic effects.aO Therefore, in most situ- 
ations it is enough to include only the term ~ hS^ in the 
Hamiltonian. 



III. BETHE-ANSATZ SOLUTION 

The most important ingredients of the algebraic 
Bethe Ansatz (BA) arc the various scattering matrices. 
The impurity-conduction electron S-matrix can be con- 
structed directly from (|]): 

i?o, = [/(f ™^(Ao - ® , (3) 

where Aq = — 1 and =^ denote the "rapidities" 
of the impurity and conduction electron j, and U is the 
[/(I) scattering matrix 

[/o,(A) = a(A)Pn + b{\)Pu + \c{X){S+ST + S+) , 

a(A) sinh(z^ + Ai?) a(A) smh{i^ + \-d) 

7{X) " sinh(i^) ' b{X) " sinh(Az9) ' ^' 
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with P|| and P-\i being the projection operators for par- 
allel and opposite spins. Excepting the small coupling 
regime, the connection of the parameters u and with 
the bare couplings Jz and J± is ambiguoustjilj and de- 
pends on the regularization procedure of the Dirac delta 
and the cutoff scheme used. Therefore, instead of Jz and 
Jj_, it is rather fj. and 9 that should be viewed as the 
basic parameters of the BA solution: fi turns out to be 
connected to the renormalized phase shift while the ra- 
tio fi/9 determines the Kondo temperature, below which 
non-Fermi-liquid correlations appear: 



TT 



(5) 



with D = N/L. 

Since electrons move with the same velocity we have 
the liberty to define their scattering matrix in a way to 
maintain integrability: 



Rij — Uij{\i — Aj) (g) Fij{Xi — Xj) 



(6) 



Here U{X) is given by Eq. (|^) and F describes scattering 
in the flavor sector: 



Xi — Xj + icXi. 
Xi — Xj + ic 



(7) 



with Xij the flavor exchange operator of particle i and j 
and c an arbitrary constant to be defined later. 

Starting from these scattering matrices we used the 
algebraic BA to determine the nested BA equations and 
then applied the dynamical fusion procedure of Ref. [ |l^] 
to eliminate the flavor degrees of freedom. The fused 
equations considerably simplify with the choice c = ^. 
Then the rapidities {Ac,; a = 1, .., M} describing the spin 
sector of the wave function satisfy 



sinh(/i(Ac 



sinli(/i(AQ 



M 

n 



sinh ^{Xa + i^) 
sinh ^{Xa — i-^) 

sinh/i(AQ, — A^ 
sinh^(AQ 



N 



X 



/3 



(8) 



The momenta of the electrons and thus the total energy 
is determined by the periodic boundary conditions 



JkALf _ 



M 

n 



sinh/x^Ao, 



N 



\ sinh^(AQ, — i^) 



E^J2f^^^ (9) 



where fkA denotes the total momentum of the fused /- 
electron composites, and L is the system size. 

In the thermodynamic limit, L, N —^ oo, N/L — 
D, the 'sp in ,, r apidities' Aq in Eq. (j^) are organized 
into stringsBlld of length r and parity v = ±: X ^ 
{Ag'''"'' ; g = 1, .., r} with 

Xi^.^) ^ xl^-"'^ = A^''^^') + [r±l ~ q] + - v) . (10) 



We have verified that to obtain a stable solution for fi < 
tt//, V andx-pmst satisfy the same stability condition as 

for / = itae 



l,..,[r/2]. (11) 



V sin(/^(7) sm{^{r — q)) > 



As shown by Takahashi and Susukijl3 the allowed (r, v) 
strings can be classified on the basis of an infinite (or 
finite) fraction expansion of (i/tt. To be specific, here 
we only discuss the simplest case ii — n/v and f < v, 
where only v different stable string configurations exist: 
n = (r, v) = (1, +), (2, -f ), .., (i. - 1, +) ani(l, -). The 
case II ^ n/f represents a singular limit :tj For / = 1 
it corresponds to the decoupling pointed while for / = 2 
it can be identified with the Emery-Kivelson point (see 
below). 



IV. THERMODYNAMICS 



A. Thermodynamic Bethe-Ansatz equations 

To derive the thermodynamic BA equations in the 
continuum limit L, N oo and D = N/L = est we 
proceeded in the usual way. We first defined the den- 
sity of rapidities (rapidity holes), gn{X) (^„(A)). These 
are related to the 'excitation energies' eri(A) through 
rin = Qn/Qn = 6^"/"^. The functions e„(A) are deter- 
mined by the following integral equations for v > f: 

e^/T ^gh v/2T ~ s * ln(l + e'-'/'^) + 6^j+ie{X) 
e^-i/T = gh v/2T + s * ln(l + e '^-^Z^) - S^j+iQ{X) 



e-j/T = s*ln 



+ ,5j,,_2 s*{l + e~'-'^) - 5j,fQ{X) {j <u~ 1). 

where s* denotes convolution with the Kernel s(A) = 
l/cosh(7rA), the driving term is given by 9(A) = 
2L>/rarctg(e'^^) and So — + —oo. The impurity contribu- 
tion to the free energy is given by 



;(A + f) ln(l + cxp(£i(A)/T)) dA , 



and, in principle, all thermodynamic quantities can be 
calculated by taking the derivatives of F""p. 



B. Analytical results 

Many of the thermodynamic properties can be deter- 
mined from the asymptotic form of the 77„'s. Using the 
Ansatz ry„(A — > ±00) ~ r]^ + b"^ ^t-^ r±\ obtains 
a set of algebraic equations for the rj^'s, 5^'s and the 
exponents t±. The latter exponents govern the scaling 
of the free energy in the vicinity of the low- and high- 
energy fixed points and are given by r+ = 4/(2 -I- /) and 



3 



T_ = 2/i/7r. The crossover between the two regimes oc- 
curs at the Kondo scale Eq. (||), which emerges naturally 
if one rewrites the thermodynamic BA equations 
in a 'universal' form by removing the cutoff Di 
The asymptotic form of the impurity free energy for 
h -^T <^ Tk is given by 



\4/2+/ 



/>2, 



^-P+(a + 6(f)0(^ln(^)) 

implying the divergence of the linear specific heat coeffi- 
cient c/T aX h — and the susceptibility as T ^ 0. The 
constants a and b above depend on the specific value of 
/ and and the-Msidual entropy S'™p is the same as in 
the isotropic casalU 



S'""P = In 



sin(^) 



sin(-^) 



(12) 



To determine the renormalization of the (/-factor and 
the Wilson ratio, we calculated the linear specific heat co- 
efficient and the bulk magnetization M*°* = —dF^°^/dh 
in the absence of the inipurity spin (but with g' — g). 
Similar to the case / = ljl3 the linear specific heat coeffi- 
cient agrees with that of the spin sector of non-interacting 
free electrons. However, the magnetization does not. Fol- 
lowing similar lines as in Ref . [ |l2| , the total spin can be 
related to ry^ and Vf^_i and thus the magnetization is 
simply given by 



d_ 



fL 



i=l 



47r 1 — ffi/n 



, (13) 



with the term, g^ fL/An, the Pauli susceptibility of free 
electrons. 

From this equation it immediately follows by integra- 
tion that at zero temperature 



F'°' = — g 



1 ,fL 



T 47rl-//i/7r 



(14) 



Thus the g-factor is renormalized due to the electron- 
electron interaction as g ^ g/{l — /i//7r)-'^/^, and to com- 
pensate the effect of Eq. (||) , we have to chose 

g ^ (1 - M//^)'/' • 

Having thus compensated the effect of the artificial 
electron-electron interaction of Eq. (H) by rescaling g, we 
can proceed to calculate the impurity contribution to the 
global susceptibility (defined with g' — \, g = 1 but no 
electron-electron interaction). We find indeed that with 
the choice g = (1— /i//7r)^/^ the low temperature (global) 
Wilson ratio, defined in terms of this global impurity 
susceptibility, takes on a universal value 



lim lim 

T->0/i 



n ^bulk imp 
^ ^glob 



(15) 



as in the isotropic case proving again that exchange 
anisottopy. is irrelevant at the two-channel Kondo fixed 
point BJ^Ea In the following we shall always denote quan- 
tities that were calculated without (i.e. compensating) 
the artificial electron-electron interaction by the super- 
script 'imp'. 

To capture the meaning of the parameter /i we also 
determined the impurity contribution to the global sus- 
ceptibility in the high temperature regime: 



imp ^ 9^ ^ 1 - /m/tt 



(16) 



Using Abelian bosonization techniques we were also able 
to prove analytically that for J± <C Jz at high tempera- 
tures 

X;Z(7^- 00)^(1-2/^)2^, (17) 

with S the phase shift generated by Jz (see Appendix ^ 
for th e der ivation and the precise definition of the phase 

shift). ElM 

This immediately implies the important relation 



II S d' 
^ = 4- -4/- 

TT n TT 



(18) 



Comparing this expression with the results of Ref. [ |2^] 
we notice immediately that /z/tt is nothing but the scal- 
ing dimension of Jj^, which satisfies the following scaling 
equation at energy scales lo well above TkS^ 



d In Jj^ ^ 

dln{ujo/Lo) TT 



(w » Tk) 



(19) 



Here luq is a high energy cutoff. For a TLS it is of the 
order of the Debye temperature, while for heavy fermions 
it is usually of the order of the Fermi energy. For a quan- 
tum dot the cutoff is the charging energy, ujq ^ Ec- The 
effective perpendicular coupling Jj_ at energy scale uj can 
be obtained by simply integrating this equation. 

Eq. (|l^) is further confirmed by noticing that for / = 2 
at the Emery-Kivelson line, S = tt/4:, the global suscep- 
tibility Eq. (^6|) identically vanishes, in complete agree- 
ment with the results of Refs. [ The point S = n/2f 
corresponding to /i — tt// is highly singular, and needs 
special care:Ej At this particular point the amplitude of 
the leading irrelevant operator, responsible for the diver- 
gence of the susceptibility, and the linear specific heat 
coefficient, becomes zero.ES 

The global susceptibility, x™ioh^ associated 

global Wilson ratio, R'^^, defined in Eq. |l^, are useful 
for magnetic Kondo systems. However, for quantum dots 
and TLS's it is rather the local impurity susceptibility 
that is of interest, i.e., the response of the system to an 
external field coupled only to the impurity spin. Hence, 
we studied the impurity contribution to the susceptibility 
when g = 1 and g' = (and in the absence of the artifi- 
cial electron-electron interaction). In order to determine 
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this we generalized the path integral derivation of Ref. [ 
[l5|] to show that 

F™P(/., T, g, g' = g) = F"^^^{h, T, g, <?' = 0) , (20) 



where h = h{l - 2f5/-K) = h{l - 
the second equality we used Eq. (| 



/7r)i/2. To establish 
Thus we can calcu- 



late the local impurity properties (those for g = 1, g' = 
and no electron-electron interaction) from the BA equa- 
tions with g' = g = {\ — f ii/tt)II? by simply rescaling 
the field h according to Eq. (pC)|).Ej 

Henceforth, unless otherwise stated, we shall denote 
the local impurity susceptibility obtained from the BA 
solution in combination with Eq. (^0|) by x""'' the 
associated local Wilson ratio by i?""P (see below). 

Following the above procedure we find at high temper- 
ature 



^imp 



1 

4T 



1 



Tk_ 
T 



(21) 
(22) 



Note that at high temperatures the specific heat exhibits 
a power law behavior which crosses over to a logarithmic 
behavior in the isotropic case /i — s- 0. Similarly, the cor- 
rections to the susceptibility about the free behavior are 
power law like and these power laws give way to logarith- 
mic corrections in the isotropic limit /i — > 0. We note that 
the above exponent_2/^/7r is formally the same as found 
for the / = 1 caseEiJ. However, the relation between /i 
and the bare couplings is quite different in the two cases 
and involves the channel number / (see Eq. (p^). 

The power-law corrections can be very easily under- 
stood from the scaling picture. Expanding the free en- 
ergy in terms of J_l one finds that the leading correc- 
tion is second order in Jj_. Making use of the scaling 
equation Eq. (|l9|) it immediately follows that the leading 
corrections to the free energy behave as T^-'^i^h ^ imply- 
ing Eq. (p^). Similar arguments lead to the conclusion 
that the impurity-induced resistivity correction behaves 
at high temperatures as 



(23) 



for T > Tk. 

At low temperatures, T <^ Tk, some care is required 
in discussing thermodynamic properties. In contrast to. 
the case / = 1 , where F is analytical around T = h = Ot2l 
here the h ~ Q, T ~ point is an essential singularity and 
the two limits T ^ and h ^ are not interchangeable. 
Taking the limit /i — s- first, we find the following non- 
Fermi liquid behavior 



X 



C" 



1 T 
T T 



(24) 
(25) 



i?""P = lim lim 



„bulk ^imp 



1 



(26) 



Here, we have defined a local Wilson ratio, i?'™P. It is ex- 
pressed in terms of the local impurity susceptibility x™p, 
and the impurity contribution to the specific heat c""p. 
It differs from the Wilson ratio defined in Eq. ^ (which 
is the usual definition for this quantity in the magnetic 
Kondo problem) by having = 0. We see that this local 
Wilson ratio depends on anisotropy /i and thus the phase 
shift (5, and is therefore not universal. 

We now consider the limit of T — > 0, with h remaining 
either finite, or taken to zero subsequently. In this case, 
and for h <^ Tk, we find from the numerical results of 
the next section, that there is a low energy scale Tfl = 

/Tk, below which the thermodynamics is Fermi liquid 
like. In particular for T <^ TpL and small magnetic fields 
h <C Tk wc find 



X™P(/i,r«TFL)'--ln(/i/TK) 

C™P(/j,r«TFL)~;^ , 



(27) 
(28) 



with TpL Ii'^/Tk (for h < Tk)- A local Wilson ratio 
for arbitrary local magnetic field h can be defined as 



i?"°P(/i) = hm 



T-,0 x'^""'(^) C™P(/l) 



(29) 



In contrast to the f — 1 case, for which this quantity 
is independent of h, but dependent, o n a nisotropy (being 
given by i?™P'f=i(/i) = 2/(1 - ^/7r)yEl), for the present 
f — 2 case it depends explicitly on both h and anisotropy 
(Fig. 1^). This important result, wiiich is consistent with 
the result for the isotropic casetHi, will be discussed in 
the following section on the numerical solution. We note 
here, however, that the local Wilson ratio for the / = 2 
case agrees with the f — 1 local Wilson ratio in the case 
of asymptotically large magnetic fields, h 3> Tk, i.e. in 
this case we have 



i?"°P(/i»TK) = 



1 



(30) 



although the meaning of /i is different in the two cases 
(see Eq. (p^)). The detailed dependence of this local 
Wilson ratio on h will be discussed in the next section. 



C. Numerical solution 

In order to obtain the thermodynamics at all temper- 
atures, it was necessary to solve the thermodynamic BA 
equations of Sec.IVA numerically. A procedure for do- 
ing this, which is valid for arbitrary values of the mag- 
netic field h and temperature has been developed in Ref. [ 
[T]| for f ~ 1. With small modifications, the same pro- 
cedure applies also to the present case. We considered 
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v=3 




FIG. 1. The impurity contribution to the entropy, 5"'"^ , 
specific heat, C™^ , and (local) susceptibility, x""''i ^ func- 
tions of T/Tk for magnetic fields H/Tk = 2"^' (solid), 
H/Tk = 2'* (dotted), H/Tk = 1 (dashed) and h/Tx = 2'^ 
(long-dashed) for the case = 3 corresponding to the largest 
anisotropy studied. 





0.8 




0.6 


H 


0.4 


s 






0.2 









0.3 


P 


0.2 


U 


0.1 







o 


1 


s 




>^ 




p 


0.5 


£ 






- 






- 
- 

h=== 




^ a 

- 




H \ \ h- 

/ \ x4 


\ \ \ 

b - 




/ \ /yy" 




— 




\ l^^-f^' 




H \ \ h- 


c 









10"' 10 * 10"' 10"' 10"' 10° lO' lO' lO' 
T/T^ 



FIG. 2. The anisotropy dependence of the entropy, specific 
heat and susceptibility at small magnetic fields {h/Tx = 2"*), 
(note that (b) is scaled by a factor 4 for comparison with the 
corresponding case of large magnetic fields shown in Fig. |^ . 
The anisotropics shown are for /^/tt = 1/v with = 3 (solid), 
V = A (dotted), u = ^ (dashed) and v = Q (long-dashed) 



anisotropies given by /i/tt = \/v with v = 3,4,5,6. In 
Fig. |l|a-c we show the thermodynamics of the anisotropic 
2CKM for a large anisotropy {v = 3) as one might have 
in a reahstic system. The characteristic non-jEermi liq- 
uid behavior, in particular the ln(2)/2 entropyllZl and the 
logarithmically divergent x™p(T) and C"°P(r)/T, are 
found at zero field. A finite field, h > 0, restores Fermi 
liquid behavior at temperatures below a low energy scale 
Tfl = h"^ /Tk, as found for the isotropic, u = oo, caset3. 
The non-Fermi liquid behavior for Q < h < Tk is there- 
fore restricted to an intermediate range of temperatures, 
Tfl < T < Tk, and we see that a clear signature of 
such behavior (such as the two peaks in the specific heat 
with each peak having only ln(2)/2 entropy, or a ln(T) 
behavior of x(T) for a temperature range below Tk) is 
possible, even at moderate magnetic fields, h ^ Tk/IQ- 

In Fig. ||a-c and Fig. ^-c we show the effect of differ- 
ent anisotropies on the thermodynamics for both a small 
external fields {h ^ Tk, Fig. g) and a large magnetic 
fields {h^TK, Fig. |). At h <^Tk, the main effect of 
anisotropy is to modify the thermodynamics at interme- 
diate, Tfl ^ T ^ Tk, and high temperatures, T > Tk- 
For large magnetic field, h 3> Tk, the anisotropy modifies 
the thermodynamic properties at temperatures, T ^ Tk 
(Fig. I). 

A characteristic feature of the present model {f — 2), is 
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FIG. 3. The anisotropy dependence of the entropy, specific 
heat and susceptibihty at large magnetic fields (h/Tx = 2^), 
(note that (c) is scaled by a factor 0.125 for comparison with 
the corresponding case for small magnetic fields shown in 
Fig. ^ ). The anisotropics shown are for fi/n = 1/v with 
u — 2> (solid), = 4 (dotted), i/ = 5 (dashed) and u = 6 
(long-dashed) 




FIG. 4. The susceptibility at high temperatures, 

T 3> Tk, for anisotropics /i/vr = 1/u and u = 3,4,5,6. 
The impurity susceptibility at T 3> Tk has the form 
^imp^y^ = ^(1/4 - B{TK/Tf^'''), i.e. the corrections to 
the free behavior 

are power laws with exponents 2^/7r = 2/3, 1/2,2/5, 1/3 for 
V = 3,4:, 5, 6. This is illustrated in the inset, which shows 
log4(rx""''(T) - 1) versus log(T/rK). Straight lines with 
slopes -2/x/7r = -1/3, -2/5, -1/2, -2/3 for i/ = 3, 4, 5, 6 are 
indicated by symbols and are well reproduced by the numer- 
ical results. 



FIG. 5. The specific heat at high temperatures, 

T » Tk, shows power law behavior C""p(T) ~ (Tk/T)^'''''' 
with fi/n — 1/p. The logarithmic derivative of this, 
dlog C'™P(T)/dlog T, (solid line) approaches ~2^i/tt at 
T ^ Tk and is shown here for v = ^ (m/'"" — 1/5)- It is 
seen to approach — 2/i/7r = —2/5 (dashed line) at high tem- 
peratures. 

that, similar to the case dissipative two state systems,0 it 
has power-law corrections in its thermodynamics at high 
temperatures (T 3> Tk)- The exponents are uniquely 
related to the anisotropy parameter ^ (cf. Eq. ^2|). This 
is in contrast to the corresponding isotropic models which 
have logarithmic corrections at high temperatures. Fig. ^ 
and Fig. || show this for the susceptibility and specific 
heat, respectively. 

Fig. H shows the remarkable magnetic field dependence 
of the local Wilson ratio, i?™P(/i). This magnetic field 
dependenceuis consistent with the result for the isotropic 
case /i O.lHI It is quite unexpected for the following rea- 
sons. First, the corresponding local Wilson ratio for the 
f — 1 case discussed in Rcfs. [ |l^ , |ll| , ^ ) is independent 
of the magnetic field and is given by 

j^imp,f=l _ 2_ — 

It depends only on the anisotropy fi (that corresponds to 
the dissipation stesngth ai in the equivalent dissipative 
two-state systemU). In contrast, for the present case 
(/ = 2), the local Wilson ratio depends both on the 
anisotropy ^ and on the magnetic field h. Thus, even 
though the regime T <^ Tfl describes a Fermi liquid, in 
the sense of Eq.^, the Fermi liquid state appears differ- 
ent to that for the / = 1 case. We note that the local 
Wilson ratio for / = 2 deviates increasingly from the 
usual / = 1 Fermi-liquid Wilson ratio with decreasing h, 
i.e., as the range over which non- Fermi liquid behavior 
dominates increases. The numerical results for i?™P(/i) 
in Fig. I indicate that i?™P(/i) vanishes as h ^ 0. This 
result is consistent with the numerical analysis of the 
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FIG. 6. The 
Wilson ratio, TC 
{y = Z (solid), V 
ai7?™P(/i) = 2 - 
set, which shows 
lines with slopes 
are indicated by 
merical results. 



magnetic field (ft) dependence of the (local) 
'^^(/i) for different anisotropies /i/vr = Xjv 
= 4 (dotted), V = h (dashed)). At /i > Tk, 
- aijKlhf''"- This is illustrated in the in- 
log(2-ff'"P(/i)) versus log(/i/TA'). Straight 
-2/i/7r = -2/3, -1/2, -2/5, for = 3, 4, 5, 
symbols and are well reproduced by the nu- 



h - 
Eq, 



estimates of the susceptibility and specific heat in 
^ In contrast, from Eq. |6| ^"°p(T ^ 0, = 0) is 
finite, so we see again that the two limits T ^ 0, = 
and T = 0, /i ^ cannot be exchanged. The inset to 
Fig. ^ shows that i?™P(/i) exhibits power law behavior 
at large magnetic fields. For h ^ Tk we find 

aii?™P = 2 - a{TKlhf^l''. 

Finally, for completeness, we show in Fig. |^ the impu- 
rity magnetization (or polarizability) in the presence of 
a local field, M^iji) — {S^){h). For each anisotropy, fi, 
Mz{h) is a universal function of H/Tk (the same holds for 
i?™P(/i)). The approach of M^(/i) to the free value for 
h ^ Tk depends on anisotropy and is found numerically 
to behave like 

M,{h) = \/2-h{TK/hf^''\ 

We expect that in the isotropic limit /i ^ the above 
power laws will give way to logarithmic corrections in the 
same way that the power law corrections to the high tem- 
perature thermodynamics gave rise to logarithmic correc- 
tions in this limit. 



V. CONCLUSIONS 

In summary, we presented a detailed analysis of the 
thermodynamics of spin-anisotropic 2-channel Kondo 
model by using the Bethe Ansatz technique combined 
with bosonization and renormalization group arguments, 
and discussed quantitatively the role of the anisotropy 
and the magnetic field. 




h/T^ 



FIG. 7. The magnetic field (h) dependence of the impu- 
rity magnetization, Mz{h,T = 0), for different anisotropies 
fi/n — {v = 3 (solid), v — 4 (dotted), i/ — 5 (dashed)). 
At h > T/f, M^{h) = 1/2 - biTK/hf"/''. This is illustrated 
in the inset, which shows log{l /2 — Mz{h)) versus log(/i/Tif). 
Straight lines with slopes — —2/3, —1/2, —2/5, for 

v — 3,4, 5, are indicated by symbols and are well reproduced 
by the numerical results. 



We showed that at high temperatures the thermo- 
dynamics is very different from that of the isotropic 
model: The local impurity susceptibility is essentially 
free impurity-like, however, the coefficient of the global 
susceptibility is non-universal, and is related to the phase 
shift 6 generated by the coupling J^. In particular, the 
global susceptibility vanishes at the Emery-Kivelson line, 
S = 7r/2/. More interestingly, the impurity specific heat 
(and the corrections to the susceptibility about the free 
behavior) exhibits a power law behavior at high temper- 
ature 



c™p(r»rK,/i) 



2^1 It: 



(31) 



The anomalous exponent /i is the anisotropy parameter 
in the Bethe Ansatz. We have shown that /x/tt is just 
the anomalous scaling exponent of the spin flip term 
at high temperature and that it is related to the phase 
shifts generated by throughEll 



4/- 



(32) 



On general scaling arguments,@@ a similar power-law 
dependence is expected to appear in the impurity resis- 
tivity. For f — 2 we find 



p™P(r » Tk) 



2/j/7r 



(33) 



rather then a simple logarithmic scaling. 

For h = 0, and for temperatures below the Kondo 
temperature, the thermodynamics is governed by the 
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isotropic 2CK fixed point and most of tlie thermody- 
namic properties resemble very much those of the fully 
isotropic model, even for strong anisotropies. 

For finite h we showed that the non-Fermi liquid be- 
havior found for ft, = persists for an intermediate region 
of temperatures TpL < T < Tk, provided h < Tk so that 
the new scale T-py^ /Tk is well below Tk, just as in 
the isotropic caseO We also showed that the Fermi liquid 
behavior below TpL is unusual in that the Wilson ratio, 
i?'"^P(/i), depended very sensitively on the magnetic field 
h. (in contrast to the / = 1 case, which is completely 
independent of h) and we calculated the detailed depen- 
dence of this quantity for several anisotropies. 

While here we focused our attention to specific values 
of the anisotropy, our calculations can be easily general- 
ized to other anisotropy values and serve as a basis for 
any interpolation necessary for cases where there is wide 
distribution of anisotropies present. 
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No. DMR 99-81283. The hospitality of the Institut-Laue 
Langevin during a visit (GZ) is gratefully . 



APPENDIX A: DERIVATION OF EQ. (17) 



To prove Eq. (|1^) let us consider the limit Jj^ ^ of 
Eq. (||) . In this case the interaction part of the Hamilto- 
nian becomes: 



i=i 



(V'|,,(0)VT,,(0)-^L,(0)Vi,,(0))5o' , (Al) 



where the electronic field operator i/)^ ^ creates chiral 
(right-moving) electron with spin a — {+, — } — {],[} 
and channel index j — {!,..,/}. Thus the interaction 
simply produces a spin-dependent potential scattering, 
and gives rise to a phase shift 5 



V'qj(x = 0+) = ^^^^(x = )e 



(A2) 



In general, the connection between Jz and 6 depends 
on the particular cutoff scheme used except for the small 
coupling limit, Jz <C 1. To be specific, here we shall use 
Abelian bosonization on a system|-af finite size L and the 
cutoff-scheme associated with itJ2J In the bosonization 
procedure we rewrite the Hamiltonian as (g ~ g' ~ 1): 



a-j j-a 



2 



2 ''^'2n 

a,j 

a 



(A3) 



where the external field h couples to the total pseudospin 
of the system, N^j denotes the total number of electrons 
with respect to the ground state with spin a in channel 
j, and the free bosonic fields satisfy: 

[5x$aj(a;),$a'j'(a;')] = 27r 6jj,Saa'S{x - x') . (A4) 
The original fermion fields can be represented as 



(A5) 



where F^j denotes the Klein factor, and a is a small 
distance cut-off of the order of the lattice spacing. 

The phase shift can be most easily calculated by intro- 
ducing charge and spin fields and quantum numbers: 



Nr. 



N, 



s,3 



(A6) 
(A7) 
(A8) 



and performing a unitary transformation on the Hamil- 
toman by U = e .j w u / v / ^ resultmg m 



Jz 

and the 'noninteracting' Hamiltonian: 



H — Hq + Hs , 



j = l fJ,=C,S 



Hs — hSn 



L 4 



-N. 



3 

N. 



(A9) 

(AlO) 
(All) 

(A12) 



From Eqs. (|A9|), (|A^), (|A5|), and (|Ag) immediately fol- 
lows that in the bosonization cut-off scheme simply 



(A13) 



To prove Eq. (|l7|) we observe that the external field 
only appears in (A12). Therefore the partition function 
factorizes as 

Z(/3) =Zo(/3) X Zs{l3,h) , 

Sg = ±l/2 

It is easy to evaluate the sum above in the L — s- 00 limit 
giving: 



g/3/ X cosh 



2 ^ ■' 4tt' 



(A14) 



The first term just generates the Pauli susceptibility of a 
free electron gas, while the second corresponds to a free 
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spin coupled to a renormalized magnetic field and gives 
a Curie susceptibility: 



^TLS 



(1 - fJj^TTf 

4T 



(A15) 



Together with Eq. ( |A13|) , this yields Eq. ([IJ). 

As a further test, one can compare this result with 
the exact relation, Eq. ([l^ ) in the small coupling limit. 
Within the Bethe Ansatz cutoff schemet3 cos(/i) = 
cos(J2/2)/ cos( J^/2), which in the appropriate small 
couphng limit gives = JJ2 + 0{Jp = 45 + 0{S'^). 
Substituting this expression into Eq. ( |l6[ ) we indeed re- 
cover the exact relation Eq. ( pT[ ) in linear order in S. 
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